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ABSTRACT 

Using gamma-ray data collected by the Astro-rivelatore Gamma ad Immagini LEggero {AGILE) satellite 
over a period of almost one year (from 2007 July to 2008 June), we searched for pulsed signals from 35 
potentially interesting radio pulsars, ordered according to °^ and for which contemporary or re- 

cent radio data were available. AGILE detected three new top-ranking nearby and Vela-like pulsars with 
good confidence both through timing and spatial analysis. Among the newcomers we find pulsars with very 
high rotational energy losses, such as the remarkable PSRB 1509-58 with a magnetic field in excess of lO'^ 
Gauss, and PSR J2229H-61 14 providing a reliable identification for the previously unidentified EGRET source 
3EG2227H-6122. Moreover, the powerful millisecond pulsar B1821-24, in the globular cluster M28, is de- 
tected during a fraction of the observations. Four other promising gamma-ray pulsar candidates, among which 
is the notable J2043H-2740 with an age in excess of 1 million years, show a possible detection in the timing 
analysis only and deserve confirmation. 

Subject headings: gamma rays: observations — pulsars: general — pulsars: individual (PSR J2229H-6114, 
PSRB1509-58,PSRB1821-24) — stars: neutron 



1. INTRODUCTION 

Although the bulk of the electromagnetic energy output of 
spin-powered pulsars is typically expected above 10 MeV, 
only ^ 0.5% of the radio pulsar population has been clearly 
identified in the gamma-ray domain (Thompson 2004; Pelliz- 
zoni et al. 2008; Halpern et al. 2008). Such meager harvest 
is to be ascribed, most probably, to the relatively low sen- 
sitivity of gamma-ray instruments with respect to radio and 
X-ray telescopes. Poor gamma-ray pulsar statistics has been a 
major difficulty in assessing the dominant mechanism which 
channels pulsar rotational energy into high-energy emission 



as well as understanding the sites where charged particles are 
accelerated. 

The dominant mechanisms and sites of the emission most 
probably depend on the rotational period P and magnetic field 
B of the neutron stars, with the millisecond pulsars (e.g. Hard- 
ing et al. 2005) behaving differently from the "classical", 
higher magnetic field ones. Population synthesis simulations, 
featuring comprehensive statistical analysis of diverse mod- 
els of emission geometry and beaming, yielded very differ- 
ent numbers of radio-loud and radio-quiet pulsars potentially 
detectable as gamma-ray emitters (see e.g., Gonthier et al. 
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2007). However, such simulations are poorly constrained by 
the source sample currently available and only a much larger 
sample of gamma-ray pulsars could help discriminating dif- 
ferent emission models. 

Astro-rivelatore Gamma ad Immagini LEggero {AGILE) is 
a scientific mission of the Italian Space Agency, dedicated 
to high-energy astrophysics (Tavani et al. 2008) launched on 
2007 April 23. The sensitivity to photons with energy in the 
range 30 MeV to 30 GeV of the AGILE Gamma-Ray Imag- 
ing Detector (GRID; Prest et al. 2003; Barbiellini et al. 2001) 
together with a time tagging accuracy of a few /is and a very 
large field of view (Ji^60° radius) all make AGILE perfectly 
suited for the detection of the time signature of new gamma- 
ray pulsars. 

The large field of view of AGILE allows long uninterrupted 
observations and simultaneous monitoring of tens of nearby 
radio pulsars belonging to the "gamma-ray pulsar region" of 
the P-P diagram characterized byfi>2xlO''G and spin- 
down energy E^i > 1 .3 x 10^^ erg s"^ (Pellizzoni et al. 2004). 

At variance with the behavior at soft X-ray energies, where 
the emission is proportional to the rotational energy loss over 
the squared distance factor {Fx °= Ed''^, where d is pulsar dis- 
tance; e.g. Possenti et al. 2002), the expected gamma-ray 
flux of radio pulsars is directly correlated to the Goldreich- 
Julian current/open field line voltage (Cheng & Ruderman 
1980; Harding 1981). It can be estimated according to the law 
oc VEd~'^ (Kanbach 2002; Arons 1996), which is known 
to reasonably fit EGRET pulsars. The large dispersion of 
such fit (probably due to different beaming fractions) pro- 
vides the minimum/maximum normalization values, allowing 
a worst/best case approach for the gamma-ray flux estimates. 
Following such an approach, we have built a sample of ~ 100 
radio-loud pulsars which are likely to be above the AGILE 
sensitivity threshold (F^in > 2 x 10"** ph cm"^ s"' at £ > 100 
MeV). Of course, we are fully aware that the actual number 
of detections will depend upon the emission model geometry 
and efficiency. 

Top-ranking targets with poor AGILE exposure, well- 
known pulsars already deeply investigated at the epoch of 
EGRET observations, and sources reserved to the AGILE 
Guest Observer Program were excluded from our list which 
encompasses 35 gamma-ray pulsar candidates.^ In this Let- 
ter we present results about isolated pulsars only. The five 
pulsars in binary systems of our sample will be discussed in 
subsequent works. 

Since gamma-ray pulsar detection must start from an at 
least approximate knowledge of recent pulsar ephemeris, a 
dedicated pulsar radio monitoring campaign has been under- 
taken. Conservatively, we present here timing analysis re- 
sults and detection claims only for pulsars having simultane- 
ous radio ephemeris for the whole relevant AGILE observing 
epochs. Our campaign will continue throughout the AGILE 
mission for most of the targets (accordingly to their visibil- 
ity), using two telescopes (Jodrell Bank and Nanfay) of the 
European Pulsar Timing Array (EPTA), as well as the Parkes 

' In ranking order: J0737-3039*, J1833-1034, J1744-1134, J174()+1000, 
J1747-2958, J2043+2740, J173()-2304, J1513-5908, J1524-5625, J1909- 
3744*, J1357-6429, J1531-5610, J1809-1917, J1617-5055, J1803-2137, 
J1801-2451, J0940-5428, J1549-4848, J1718-3825, J1824-2452, J1730- 
3350, J0900-3144*, J142()-6048, J1739-3023, J1751-2857*, J1804^2717*, 
J2229+6114, Jl 105-6107, J1721-2457, Jl 124-5916, J1722-3712, J1740- 
3015, J1823-3106, J1745-3040, J1016-5857 (the asterisk indicates binary 
systems). See http://agile.asdc.asi.it for details about AGILE Data Policy and 
AGILE Team Targets List. 



radio telescope of The Australia Telescope National Facility 
(ATNF). Details about radio data analysis procedure in use for 
AGILE are described in Pellizzoni et al. (2008). Timing solu- 
tions for the detected pulsars can be requested by contacting 
the corresponding author. 

2. AGILE OBSERVATIONS AND DATA ANALYSIS 

Pulsar data were collected since the early phases of the mis- 
sion. Timing observations suitable for pulsed signal analysis 
started in 2007 July (at orbit 1 146) after engineering tests on 
the payload. In this Letter, we analyze data collected up to 
2008 June 30. 

AGILE pointings consist of long exposures, typically last- 
ing 10-30 days, slightly drifting (<1° day"^) from the start- 
ing direction in order to match solar panel illumination con- 
straints. The AGILE total exposure (computed with the GRID 
scientific analysis task AG_ExpmapGen) for the 11 month 
data span considered is typically <10^ cm^ s for each target. 

Data screening, particle background filtering and event di- 
rection and energy reconstruction were performed by the AG- 
ILE Standard Analysis PipeHne (BUlLD-15). We adopted the 
AGILE event extraction criteria and timing procedures cali- 
brated and optimized with the observations of known garmna- 
ray pulsars as described in Pellizzoni et al. (2008). In par- 
ticular, we performed our timing analysis looking for pulsed 
signals using both £ > 50 MeV G class events (i.e. events 
identified with good confidence as photons) and the combina- 
tion of Gh-L events (L events are significantly contaminated 
by particle background) collected within 60° from the center 
of GRID field of view and with an extraction radius around 
pulsar positions optimized as a function of photon energy ac- 
cording to the point-spread function (PSF) and background 
level of AGILE. For noisy observations in confused regions 
(e.g. for J2229-H6114), an extraction radius of ~2° (well be- 
low £■ > 50 MeV PSF^) was used to optimize the signal-to- 
noise ratio. Source position accuracy is <0?5 at present level 
of targets exposure. 

Given the measured AGILE'S time tagging accuracy of 
<200 /Its and the good radio monitoring (i.e., valid epoch 
range, adequate number of time of arrivals (ToAs)) available 
for the majority of our targets, the most significant pulsed sig- 
nal detection is typically expected within the errors of the ra- 
dio ephemeris frequency values. In particular, we performed 
standard epoch folding^ over a frequency range defined ac- 
cording to 3(7 errors of radio ephemeris. 

The frequency resolution is 0.5/Tjs (where 7"^^ is the time 
span relative to the coverage of each target) and we oversam- 
pled it by a factor of 10 {0.05/Tjs is the frequency step of our 
grid). Since the AGILE useful data span for pulsar detection is 
typically of the order of few months, implying a period search 
resolution of <10~"^ s, the resulting number of independent 
period search trials turns out to be few tens while no indepen- 
dent trials on period derivatives are needed. 

Pearson's x statistics is applied to the 10-bin folded pulse 
profiles resulting from each set of spin parameters, yielding 
the probabilities (weighted for the number of trials performed 
on the data set) of sampling a uniform distribution, assessing 

2 The on-axis AGILE PSF 67% containment radius is ^ 10° for £ = 50 
MeV, ~ 5° for E = 100 MeV, and 0?5 for £ = 1 GeV. 

' When adequate radio observations covering the time span of the gamma- 
ray observations are available (i.e. where the WAVE terms are included in 
TEMP02 ephemeris files; (Hobbs et al. 2004, 2006; Edwards et al. 2006)), 
we account also for the pulsar timing noise in the folding procedure as re- 
ported in Pellizzoni et al. (2008). 
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TABLE 1 

Emission parameters of the seven pulsars discussed in text (Firmly detected gamma-ray pulsars are reported in the upper part 

OF the table) 



PSRName G.Lon. G.Lat. P IT \ogE xldi^stf Tume" (^space^ ^V" logL^^ L^/E 

deg deg ms yr kpc erg s~' erg s"' 

J2229+6114 106.65 2.95 51.6 1.0x 10"' 12.0 37.35 6.0(36) 5.0 7.5 26±4 35.36 0.01 

J1513-5908 320.32 -1.16 150.7 1.6x 10^ 5.8 37.25 4.2(3) 4.0 6.4 34±6'' 35.04 0.006 

J1824-2452 7.80 -5.58 3.0 3.0x 10^ 4.9 36.35 4.2(1) 4.2 3.6 18 ±5 34.62 0.02 

J1016-5857 284.08 T88 107.4 2.1 x 10" 93 36.41 6.0(69) 48 123 62 ±6* 35.71 (12 

J1357-6429 309.92 -2.51 166.1 7.3 x lO' 4.1 36.49 5.2(7) 4.7 1.8 <14 <34.35 <0.007 

J2043+2740 70.61 -9.15 96.1 1.2x 10* 1.1 34.75 4.1(1) 4.2 0.6 <6 <32.84 <0.01 

J1524-5625 323.00 0.35 78.2 3.2x 10" 3.8 36.51 4.6(4) 4.3 1.0 <16 <34.34 <0.007 



^Dispersion measure distance and characteristic age obtained from ATNF Pulsar Catalog (http://www.atnf.csiro.au/research/pulsar/psrcat/; Manchester et al. 
2005). For J2229+6114, Halpem et al. (2001) derived a much shorter alternative distance of 3 kpc from X-ray absorption. 
''Pearson's statistics appUed to the 10-bin, £ > 50 MeV folded pulse profiles (iVj, is the number of steps over period search grid, see text). 
"^Pulsed emission detection significance (weighted accounting forNst). 

"•Likelihood analysis detection significance and flux (or 2a upper limit) in units of 10"* ph cm""^ s"' (£ > 100 MeV), obtainined by search for a source near the 
pulsar direction. 

^ £ > 100 MeV luminosity assuming traditional 1 sr beam emission and photon index 2.0 (for J2229+61 14 the measured photon index of 2.2 has been used). 
'Possibly multiple sources. Source confusion may affect flux and significance. 



the significance of the pulsed signal (sinusoidal or not). In 
fact, applying this statistics to the well-known EGRET pulsars 
(PeUizzoni et al. 2008), we obtained a perfect match between 
the best period resulting from AGILE gamma-ray data and 
the period predicted by the radio ephemeris with discrepan- 
cies comparable to the period search resolution. Furthermore, 
we verified our timing results also applying bin-independent 
parameter-free statistics as the test (Buccheri et al. 1983) 
and the //-test (de Jager et al. 1989) that are typically more 
sensitive than tests for the search of sinusoidal pulses. In 
particular, the use of the //-test is suggested if no a priori in- 
formation about the light curve shape is available as in our 
case. 

First, for each target we searched for pulsed signals using 
the whole available data span. Later, each observation block 
was analyzed to check for possible flux and/or pulse profile 
variability. 

We also performed a prehminary maximum likelihood 
analysis (ALIKE task) on the AGILE data for the regions con- 
taining our targets in order to exploit the instrument's imag- 
ing capabilities to assess gamma-ray source parameters. Here, 
we focus on timing analysis leaving detailed source position- 
ing, flux and (phase-resolved) spectral analysis for future pa- 
pers when higher counts photon statistics on each target will 
be available and instrument effective area caUbrations will be 
consolidated. 

3. NEW GAMMA-RAY PULSARS 

Table 1 Usts the emission parameters of the seven pulsars 
discussed in this Letter The resulting radio-aligned light 
curves are plotted in Figure 1 where for each pulsar the ac- 
tual bin size and energy range has been adjusted according 
to the available statistics and hght-curve structure. We note 
that in all cases radio and gamma-ray timing results are com- 
patible, with the highest significance frequency detected in 
gamma rays within the errors of the radio ephemeris value, 
considering also the period search resolution. Examples of 
the exploration of much larger frequency search grids are 
shown in Figures 2 and 3 for PSRJ2229-h6114, for which 
the best gamma signal is within Icr from the radio peak, and 
for PSR J 1824-2452, for which the radio-gamma frequency 
discrepancy is comparable with the gamma-ray period search 



resolution. Both statistics and Z^///-test provide compara- 
ble detection significances, except for PSR J2229H-6114. The 
test appUed to this pulsar provides slighly better results 
than the other statistics (Z^ = 39.5, Zj = 45.9 correspond- 
ing to an ^6a detection). Furthermore, we verified that our 
analysis procedure (potentially affected by instrument-related 
systematic errors and biases in events extraction criteria) does 
not produce fake detections at a significance level above 3a 
when the radio-ephemeris are applied to randomly extracted 
AGILE data. 

Four targets are also firmly detected by the likelihood spa- 
tial analysis. It is not surprising that spatial detection be miss- 
ing for weaker targets at the present AGILE exposure level. 
In high-background regions of the Galactic plane, our pulsed 
emission search sensitivity can be better than our spatial anal- 
ysis sensitivity. Furthermore, timing analysis can be applied 
to the full AGILE event list (i.e. Gh-L event classes extending 
up to 60° from the center of the field of view), while, cur- 
rently, spatial analysis is only fully calibrated for the G class 
events detected within 40° from the center of the field of view. 

However, some of the non-detections in spatial analysis are, 
indeed, puzzling. According to Figure 1, PSR J2043H-2740 
has a pulsed flux and Ught curve very similar to PSR J 1513- 
5908 (note that both are based on Gh-L events). Moreover, 
J2043H-2740 is out of the Galactic plane, in a region well ex- 
posed and with low diffuse emission. Thus, it should have 
been detected by the image analysis more easily than J1513- 
5908, which sits in a higher background region, unless, of 
course, they have very dif^ferent spectra. 

Pulsars in Table 1 are ranked according to their overall de- 
tection significance. Pulsars firmly detected both through tim- 
ing and spatial analyses are placed above the line. Pulsars 
detected through timing analysis alone await longer simulta- 
neous gamma-ray and radio observations for confirmation and 
reliable luminosity estimation. 

In Table 1, the most significant detection is 
PSR J2229H-61 14 for which AGILE detected pulsed emission 
(radio/gamma-ray periods discrepancy <10~'^ s, Figure 
2) and pinpointed the most likely position (/ = 106?86, 
b = 2.°94) to ^0.°2 from the radio pulsar. Our detection pro- 
vides a reliable identification for the previously unidentified 
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Fig. 1. — Background-subtracted E > 100 MeV folded pulse profiles (ex- 
cept for the E > 50 MeV J 1357-6429 light curve) of the pulsars shown 
in Table 1. Events extraction radius is optimized as a function of energy 
(see the text). All pulse profiles are obtained using G event class only 
except for PSRJ15 13+5908 and PSR J2043+2740 that have higher counts 
statistics because detections are obtained including G+L event class. For 
PSRJ2229-I-6114 the histogram of the --75 (G+L) events with E > 1 GeV is 
also shown. Absolute timing is performed for each target: the main radio 
peak (1.4 GHz) corresponds to phase 0. Possible fluctuations of the dis- 
persion measure over the considered time interval are not expected to sig- 
nificanfly affect phasing results, given the time resolution of the available 
gamma-ray light curves. 



EGRET source 3EG2227H-6122 (/ = 106.°53, b = 3.°18). 
The AGILE source position, pulsed flux, and photon index 
(~2.2) are consistent with the EGRET values (Hartman et al. 
1999). The gamma-ray light curve of this pulsar (detected up 
to over 1 GeV), featuring just one prominent peak shifted ^ 
180° in phase from the radio main peak, is shown in Figure 
1. It is worth noting that assuming a distance of 3 kpc 





-6x10-" -4x10-" -2x10-" 2x10-" 4x10-" 6x10-" 8x10-" 
P^GILE-P_RADIO [s] 

Fig. 2. — PSRJ2229+6114: gamma-ray period search result (period 
trials versus Pearson statistics). The radio period (vertical line at 
i'AGiLE--f^iADio = 0) is 51.64101208(3) ms (PEPOCH = 54575.5856 MJD). 
The most significant gamma-ray pulse profile is obtained for Pf^,',p = 
51.641012085(2), a value within the la error on the radio period (vertical 
dashed Unes). 



inferred by X-ray observations (Halpem et al. 2001) and 
isotropic emission, the pulsar gamma-ray efficiency would 
be ~ 0.5, a factor of 20 higher than that quoted in Table 1. 
In modern pulsar beaming models (in particular high-altitude 
models) the assumption of isotropy in luminosity calculation 
could be a better approximation in most cases, implying 
that efficiencies should be increased by roughly an order of 
magnitude from those based on the commonly assumed 1 sr 
beam (Walters et al. 2008). Anyway, precise efficiency mea- 
surements await for a better assessment of pulsar distance, 
flux, and beam geometry. 

PSRJ1513-5908 (B1509-58) was detected by COMPTEL 
in the 1-10 MeV range, while EGRET reported (Fierro 1996) 
only marginal evidence for a weak <4(t source at ^ 1° from 
the radio position, with a pulsed emission upper limit of 
<58 X 10 ph cm~^ s~'. The AGILE discovery of pulsed 
emission from PSRB 1509-58 relies also on a 6.4cr detec- 
tion of a gamma-ray source (possibly multiple sources) at 
<0?4 from the radio position. The main gamma-ray peak 
at phase ^ 0.35 in the AGILE light curve is aligned with the 
soft gamma-ray peak seen by COMPTEL in the 0.75-30 MeV 
band (0.35 ± 0.03; Kuiper et al. 1999) shghtly traihng flie hard 
X-rays single peaked light curves (Saito et al. 1997; Rots et al. 
1998). A second possibile peak in the AGILE light-curve at 
~ 0.85 could have correspondence with the marginal feature 
seen in the 10-30 MeV COMPTEL fight curve. The COMP- 
TEL sharp spectral break between 10 and 30 MeV (Kuiper 
et al. 1999) is confirmed by our AGILE point at 100 MeV. 
Our data imply a softening of the photon index from ^1.7 
to ^ 2.5 going from tens to hundreds of MeV. Such a low- 
energy break, compared to the more common GeV spectral 
break in gamma-ray pulsars, could be the telltale signature of 
the photon splitting process inducing significant e.m. cascade 
attenuation due to the strong magnetic polar field (greater than 
10'^ G) of this pulsar (e.g. Harding et al. 1997). 

At variance with all the other targets, the millisecond pul- 
sar J 1824-2452 in the Globular Cluster M28 was detected by 
AGILE with good significance (greater than 4a) and perfect 
radio-gamma periods match, only in the time interval 54339- 
54344 MJD (Figure 3). The main radio peak at 1.4 GHz is 
coincident with the broad single peak seen in gamma rays. 
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-1x10-1" -5xlO-'l 5xlO-'l 1x10-1" 

PJ\GILE-P_RADIO [s] 

Fig. 3. — PSRJ1824-2452: gamma-ray period search result. The radio 
period is 3.0543151208713(1) ms (PEPOCH = 51468 MJD). The most sig- 
nificant gamma-ray pulse profile is obtained forP^g^LE = 3.05431511(1) ms. 
The corresponding radio-gamma periods discrepancy <10~" s is compara- 
ble to the period search resolution in the considered data span (54339-54344 
MJD). 

Only marginal detection was obtained integrating other ob- 
servations with comparable exposure or the whole data span. 
Noise fluctuations could possibly explain the apparent vari- 
abiUty. Alternatively, although its gamma-ray efficiency and 
high stability of spin parameters are compatible with rotation- 
powered emission, some additional mechanism disturbing the 
neutron star magnetosphere in the dense cluster environment 
could be invoked to explain the variable gamma-ray phe- 
nomenology of this peculiar pulsar. AGILE timing failures at 
submillisecond level in some observations (mimicking source 
variability) cannot be excluded. However, this seems unlikely, 
since we verified both timing accuracy and stability at ^ 200 
IIS level with Vela pulsar observations. Confirmation of this 
tantaUzing result about physical variabiUty will rest on longer 
monitoring campaigns. 

PSR J1016-5857 stands out in Table 1 for its very high effi- 
ciency in converting rotational energy loss into high-energy 
gamma rays. This may be ascribed to distance uncertain- 
ties, as happened to be the case for the recently discovered 
PSR J2021-H3651 (Halpem et al. 2008) for which the dis- 
tance derived from the dispersion measure is certainly overes- 
timated. Furthermore, we note that the position (I = 284P47, 
b = -0?94) and flux of the AGILE gamma-ray source is only 



marginally compatible with 3EG 1013-5915, for which multi- 
ple associations were proposed (Hartman et al. 1999). The re- 
gion, originally covered by the gamma-ray source 2CG 284- 
00, discovered by COS-B (Swanenburg et al. 1981), is very 
complex. Other nearby sources, and in particular pulsar 
PSR J1016-5819 (not belonging to our sample), could signif- 
icantly contribute to the high flux observed by AGILE in this 
region. 

Of the remaining three pulsars, detected only through tim- 
ing analysis, the most notable is certainly PSR J2043-I-2740 
which, with an age in excess of one million years, would be 
the oldest ordinary (nonrecycled) pulsar seen in the gamma- 
ray domain. The precise measurements of its flux will yield its 
luminosity and thus its efficiency, a parameter of paramount 
importance for the understanding of ganmia-ray emission as 
a function of the pulsar age. 

According to their timing and spatial analysis outcomes, 
Ught-curve profiles, and radio-gamma phasing upshots, 
PSRJ2229H-6114, J1513-5908 and J1824-2452 can be con- 
sidered rather solid detections. If the detection of J 1016- 
5857, J2043-I-2740, J1357-6429, and J1524-5625 will also be 
firmly established, our seven object sample, together with the 
detection of J2021H-3651 (Halpern et al. 2008), would imply 
that gamma-ray emission is a common feature of luminous 
radio-loud pulsars, be they young or old. Indeed, our list en- 
compasses the second youngest (J1513-5908) and by far the 
oldest nonrecycled pulsar (J2043H-2740) detected at gamma- 
ray energies. The gamma-ray pulsar search and the improve- 
ment of AGILE timing procedures requiring unprecedented 
long integrations is an ongoing effort, in close collaboration 
with the radio observers. We trust that the new gamma-ray 
pulsar detections presented here will trigger multiwavelength 
observations (e.g. in X-rays) simultaneous to the continuously 
improving AG/LE and Fermi (formerly GIAST) exposures on 
these targets. The availability of high-resolution light curves 
and their precise multiwavelength phasing wiU pave the way 
to improve pulsar models as well as our understanding of the 
pulsar population as a whole. 

We thank the referee for her/his useful recommendations 
and comments. The Parkes radio telescope is part of the 
Austraha Telescope, which is funded by the Commonwealth 
of Australia for operations as a National Facility managed 
by CSIRO. We thank Jules Halpern and Fernando Camilo 
for useful discussions and contributions on gamma-ray pul- 
sar timing. 
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